Observation of local temporal correlations in trapped quantum gases 
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We measure the temporal pair correlation function (r) of a trapped gas of bosons above and 
below the critical temperature for Bose-Einstein condensation. The measurement is performed in 
situ using a local, time-resolved single-atom sensitive probing technique. Third and fourth order 
correlation functions are also extracted. We develop a theoretical model and compare it with our 
experimental data, finding good quantitative agreement and highlighting the role of interactions. 
Our results promote temporal correlations as new observables to study the dynamics of ultracold 
quantum gases. 
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The intriguing effect of particle bunching was first ob- 
served in a seminal experiment by Hanbury Brown and 
Twiss (HBT), where they studied correlations among 
pairs of photons coming from a chaotic source [l| . The re- 
sult, which could be explained in terms of classical waves, 
had a difficult route to be accepted under a particle per- 
spective. The full quantum theory, due to Glauber 
signed the birth of quantum optics and made the formal- 
ism, with all its physical content, available for massive 
particles. Over the years, analogous HBT experiments 
were performed with electrons [3|, neutrons [4[ and cold 
atoms Q . The possibility to extract information about 
the quantum statistics and the coherence, paved by the 
HBT experience, conjugated with the capabilities of de- 
riving the temperature and the spatial order [7[, con- 
tributed to make this technique one of the most powerful 
to probe atomic systems. The quest to understand the 
behavior of more and more complex samples suggests its 
application to the study of strongly interacting ID gases 
disordered [10] and supersolid phases jll| and to 
identify non trivial excitations [12|- While first order 
correlations are often accessible via interference exper- 
iments, higher order correlations require in general the 
recording of density or atom number fluctuations by a 
probe sensitive enough to detect single particles (count- 
ing techniques) or, at least, atomic shot-noise (absorption 
imaging). In order to have a good statistical description, 
an average over many realizations of the system (in the- 
ory all possible realizations) is needed. Consequently cor- 
relations, especially at orders higher than two, are usually 
difficult to measure because of the huge statistics required 
for a reliable signal. Only in some limited cases, intrinsic 
processes in a quantum gas such as photoassociation or 
three body losses can be used as a sensitive probe for 
higher order correlations at zero distance [HI, Q . The 
direct observation of third order correlations is still chal- 
lenging and, using standard techniques, requires a con- 
siderable effort in data collection and analysis [T^j. In 
this, like in the great majority of the above mentioned 
experiments with ultracold gases, correlations have only 
been studied in the spatial domain whereas the temporal 
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FIG. 1. (Color online) a) Schematic of the apparatus. The 
electrons of the focussed beam collide with atoms leading to 
ionization with a probability of 40%. The ions are guided 
to the detector to be counted. The average count rate is 1.7 
kHz. b) The electron- atom interaction volume is defined by 
the electron beam itself and by the extension of the cloud 
along the z axis, c) A typical measurement of the ion signal 
showing the first 20 bins. 



counterpart has been very poorly explored, limited only 
to the characterization of atomic beams [5|, [l6| . Boosted 
by recent achievements [TtI-[iq| . spatially resolved single 
atom sensitive detection methods can give direct access 
to higher order correlation functions in situ. Such experi- 
ments are not complicated by the time of flight expansion 
and directly manifest the properties of the many-body 
system. Moreover if the detection method is only locally 
probing the system, even time-resolved correlation mea- 
surements are possible. 

In this Letter, we report on the first observation of 
temporal thermal bunching of ultracold bosonic atoms 
in a trap. Using scanning electron microscopy as time- 
resolved local detection method, our measurements di- 
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rectly yield the second, third and fourth order time cor- 
relation functions. Notably, due to the spatially resolved 
measurement and due to the possibility to follow the dy- 
namics of the system, our technique is effectively an orig- 
inal method for characterizing ultracold quantum sys- 
tems, both in space and time. All these features can turn 
out to be extremely useful in providing deep insight into 
the dynamics of strongly correlated many-body quantum 
systems [2o| . 

The general form of the normalized spatio-temporal 
correlation function of n particles at position at time 
ti, with i = 1, n, is given by: 



(n) (ri,*i;...;r n ,t n ) = 



(1) 



where \I/ are the bosonic operators and (...) indicates the 
ensemble average. We first derive an analytical expres- 
sion of g( 1 \ri,ti]T2,t2) for an ideal Bose gas at tem- 
perature T, trapped in a harmonic potential V(r) = 
mcj 2 r 2 /2, with average trapping frequency lj, extending 
the approach of Ref. [21] to take into account also the 
temporal evolution. Given r = — t\ and r = Y2 — 1*1 
and assuming cjt, hu/(kBT) <C 1, we obtain: 
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where r c = is defined as the correlation time. From 

kbj- 

the above expression we can derive any higher order cor- 
relation function for thermal bosons and, in particular, 
the second order correlation can be easily calculated as 
g( 2 \r, r) = l + |^^^(r, r)| 2 . The time correlation function 
(0, r) can be interpreted as the probability to detect a 
particle a time r after another particle at the same posi- 
tion (r = 0). For a thermal cloud of bosons this function 
decreases from 2 to 1 on a time-scale related to the cor- 
relation time r c . A value of the pair correlation function 
higher than 1 indicates bunching of thermal bosons in 
time. For a coherent source like a BEC, we expect in- 
stead # (2) (0,r) = 1 for any r [2l|, i.e. a flat detection 
probability in r. 

Since our setup has been described earlier 22, 23|, we 
briefly illustrate the experimental procedure we have fol- 
lowed to measure the atomic correlation functions. The 
key feature of our experiment is a scanning electron mi- 
croscope which is implemented on a standard apparatus 
for the production of ultracold quantum gases. The mea- 
surement principle is based on the electron impact ioniza- 
tion of the atoms with subsequent ion detection. A sketch 
of the working principle is depicted in FigHJ In addition 
to a spatial resolution of better than 150 nm, the tech- 
nique is characterized by a sequential detection method. 
Thus, time-dependent quantities such as the second or 
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FIG. 2. (Color online) Normalized second order temporal cor- 
relation function. Data (dots) acquired at different tempera- 
tures (circles for 45 nK, squares for 100 nK) are plotted to- 
gether with the fitting functions (lines) explained in the text. 
The inset displays the data (diamonds) acquired for a BEC 
and the corresponding fit (line). Please note that even well 
below the critical temperature (the thermal fraction cannot 
be detected in time of flight absorption imaging) we are able 
to measure a small residual bunching induced by the thermal 
component. 



higher order correlation function g^ n \r), n = 2,3, .., be- 
come experimentally accessible. To prepare the atomic 
sample we load approximately 1.5 x 10 6 87 Rb atoms from 
a magneto-optical trap in a CO2 dipole trap. The atoms 
are then evaporatively cooled above or below the criti- 
cal temperature. We can typically prepare cold thermal 
samples of 2 — 3 x 10 5 atoms at 150 nK as well as BECs of 
up to 1 x 10 5 atoms in the F = 1 hyperfine ground state. 
Since integration along the probing line reduces the cor- 
relation signal amplitude (see for example [6]), after the 
evaporation we compress the cloud along the direction of 
the electron beam by adiabatically transfering it to a light 
sheet dipole trap. This trap is realized by means of a fo- 
cussed elliptical 852 nm laser beam, with waists (6.5, 130) 
fim. We then completely switch off the CO2 trap and the 
sample is held by the light sheet alone for all the duration 
of the measurements (Fig. [1]). The final frequencies, for a 
beam power of 5 mW, are oj XiViZ ~ 2tt x (13, 27, 580) Hz. 
To detect the atoms, an electron beam of 6 keV energy, 
20 nA current and 120 nm FWHM is focussed at the cen- 
tre of the cloud, where the atom density and hence the 
number of ions produced is maximal. This choice is only 
meant to allow for a reduced number of experimental cy- 
cles but any other point of the cloud can be investigated 
as well. The electron-atom interaction is confined within 
a distance smaller (x, y directions) or comparable (z di- 
rection) to the typical correlation lengths of our atomic 
samples. The ions produced by electron impact ioniza- 
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tion are collected in a channeltron. For each bin of 10 /is 
we record the number of ions detected and the absolute 
detection time, as shown in the inset of FigJTJ In a typical 
measurement of 600 ms duration we extract ~ 1000 ions. 
After the measurement, the remaining atoms are released 
from the trap and imaged after time-of-flight (TOF) to 
extract the temperature. In order to optimize signal- 
to-noise ratio we compute the correlation functions over 
about 1500 repetitions of the experiment. 

From Eq.(pQ), integrating on the absolute time and av- 
eraging on the different repetitions of the experiment, 
we calculate the correlation function g^ 2 \r) for a cold 
thermal cloud. In Fig|2]we report the measurements for 
two different temperatures T = 45 nK and T = 100 nK. 
In both cases, we can clearly observe bunching. In these 
data plots we have omitted the first two bins, which show 
an abnormally high value of g^ 2 \r). These points are af- 
fected by extra ion counts introduced by the electronics 
of the detector. As a test for the procedure, we apply 
the same technique to a reference measurement which is 
a priori uncorrelated. This is obtained by probing a ther- 
mal cloud at T = 230 nK with an electron beam of waist 
well larger than the correlation length. As expected, no 
bunching signal is detectable in the measurement, mak- 
ing sure that what we have observed in the data sets 
taken in standard conditions is genuine thermal bunch- 
ing. 

Averaging over several experimental cycles, fluctua- 
tions in the total atom number affect the normalization 
of the correlation function. As a result, an offset shifts 
the uncorrelated signal to a value 1% above 1. To com- 
pensate for these fluctuations we normalize g^ 2 \r) by 
the factor 1 + (a 2 - (N))/(N) 2 , where a 2 and (N) are 
respectively the variance and the mean value of the total 
atom number in the different experimental realizations. 
A measurement with a BEC is also presented in the in- 
set of FigJ21 Notably, we can still detect a small corre- 
lation signal due to a residual thermal fraction present 
well below the critical temperature T c . The g^ 2 \r) data 
points are fitted with the function we derived for an ideal 
non interacting gas of bosons integrated over the volume, 
leaving the amplitude and r c as free parameters (lines in 
Fig|2j). The results of the fits are shown in Fig|3j where 
the extracted correlation times are plotted as a function 
of the temperature of the sample. Since we cannot de- 
rive the temperature for a cloud with more than 70% 
condensate fraction by TOF measurements, in the plots 
of Fig|3] we indicate as the BEC temperature the aver- 
age temperature between zero and the coolest measurable 
temperature of 30 nK. The agreement between the exper- 
imental data and the theoretical function r c = h/iksT) 
is fairly good. Hence this technique can be proposed as 
a local probe for the temperature of ultracold samples, 
especially when standard imaging techniques fail. 

In the inset of Fig|3]we show the fitted amplitude of the 
normalized second order correlation functions at different 
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FIG. 3. (Color online) Correlation time r c as a function of the 
temperature. The values obtained fitting the experimental 
data (circles) are compared with the function r c = h/(kBT) 
(line). The inset shows the fitted bunching enhancements 
(circles) and the corresponding theoretically expected values 
for the non-interacting (solid line) and the interacting (dashed 
line) model. 




FIG. 4. (Color online) Normalized third order correlation 
function (dots) along the axis (n = £2 — £i,T2 = £3 — £1 = 2ri) 
for a thermal cloud with T — 100 nK. The solid line is the 
fitting function explained in the text. 

temperatures together with the values expected from the 
volume integration of the non-interacting model and of 
its extension to the interacting case [2l| : 

^(r,0) = l + ^ + | 5 W(r,0)| 2 (l-^,) (3) 

being a the s-wave scattering length. As expected [2l|, re- 
pulsive interactions play a significant role for in situ mea- 
surements since they induce short-range anti-bunching 
that reduces or overcomes the bunching signal of bosons. 
Our measurements are well reproduced by the interacting 
model. In the regime where the BEC fraction is domi- 
nant, the signal is further reduced by the presence of a 
large coherent component in the integrating volume and 
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by the role of interactions in the condensate. 

In addition, we also calculate the third order correla- 
tion function, extracting it from the same data sets for 
which we derived the g( 2 \r). Bunching is expected to 
be more pronounced at higher orders n as a consequence 
of the factorial law n\ that regulates the dependence of 
the correlated to the uncorrelated amplitudes. For this 
reason higher order correlations can be employed as a 
highly sensitive test for coherence, with the only draw- 
back represented by the need of high statistics, as proved 
by the very few experiments reporting correlations at or- 
ders higher than 2. 

In Fig|4]we show an example of the third order corre- 
lation function ^^ 3 ^(ri,T2), obtained for the data set at 
T = 100 nK, along the axis (n, T2 = 2ti). The solid line 
in the picture is a fit made with the volume integration of 
the non interacting model g^\r\,T2) = 1 + |<7^( r i)| 2 + 
\9 {1 Hr2)\ 2 + \9 {1 Hr2-T 1 )\ 2 +2^gW(T 1 )g^(-T 2 )g^(T 2 - 
ti)), along the same axis, leaving the amplitude as the 
only free parameter. The measured amplitudes of the 
third order correlation functions g 3 (0, 0) — 1 are 0.14±0.07 
and 0.10 ± 0.02, respectively for the data sets at T = 45 
nK and T = 100 nK. Our measurements show that the 
amplitudes of the correlations are also affected by 
repulsive interactions since the bare volume integration 
of the non-interacting model is not sufficient to explain 
the observed reduction of the signal. Interestingly, the 
corresponding ratios (g 3 (0,0) — l)/(g 2 (0) — 1), which are 
3.2 ± 2.4 and 5.0 ± 2.2, show instead agreement with 
the values derived from the non-interacting model: 4.2 
and 2.8. This result is in accordance with the semi- 



nal work on the three-body losses [14[ and may signal 
that the effect of interactions does not scale with the or- 
der of the correlation function. Finally we measured the 
fourth order correlation function amplitude: for T = 45 
nK, g^ (0,0,0) - 1 = 0.46 ± 0.42 and for T = 100 nK, 
# (4) (0,0,0) - 1 = 0.23 ± 0.05. At times longer than r c , 
the g^ correlation signals drop respectively to the val- 
ues of 0.8 ± 0.1 and 1.03 ± 0.01. The evaluation of g^ 
correlations turns out to be strongly affected by the poor 
statistics, thus rendering any comparison with theory dif- 
ficult. Further investigation on fourth and higher order 
correlations can be, however, useful to identify the influ- 
ence of interactions in multiple-particle bunching. 

In summary, we have observed temporal pair correla- 
tions in a cold gas above T c . We have measured g^ 2 \r) 
for different temperatures of the atoms and compared 
the results with a theoretical model, that we derived for 
a non-interacting and for an interacting system of ther- 
mal bosons. A measurement on a BEC has also revealed 
a minimal bunching compatible with the presence of a 
residual thermal fraction. A signal of third and fourth 
order correlations has been derived from the same data. 
The amplitudes of the bunching enhancement of g^ and 
g( 2 ) correlations are affected by repulsive interparticle in- 
teractions, while their ratios are compatible with a non- 



interacting modellization. The electron microscopy tech- 
nique, which we use, represents a powerful, fast and effi- 
cient way to measure correlations, with the possibility to 
probe the system in space and time. Notably, pair corre- 
lations in time can give access to the dynamical structure 
factor and reveal, in contrast to static pair correlations, 
the dynamical properties of the quantum system. 
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